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I. INTRODUCTION 



The combined effect of carrier interactions and Pauli 
principle leads to itinerant ferromagnetism in metals 
whenever the kinetic energy gain by parallel alignment 
of electronic spins is exceeded by the exchange energy 
of the antisymmetric wave function ( lj). The simplest 
model is the Hund's rule in atoms, where the energy is 
minimized when the fermionic wave function corresponds 
to the alignment of spins. Meanwhile the electron-spin 
correlations seems to be an important factor in the two- 
dimensional (2D) metal-insulator transition :2j. Mea- 
surements of the in-plane magnetoconductivity of a dilute 
2D electron system in silicon heterostructures have shown 
the evidence for a zero-temperature quantum phase tran- 
sition at critical density, indicating the existence of a fer- 
romagnetic instability (J3,H)- On the other hand a new 
method that measures the minute thermodynamic spin 
magnetization of a dilute 2D electron system in a silicon 
inversion layer favors the paramagnetic phase over the 
spontaneous magnetization 

The interaction strength between electrons is tradi- 
tionally described by the 2D Wigner-Seitz radius, r s = 
l/(7rn s ag) 1 / 2 , where as is the Bohr radius and n s = 
n/L 2 sheet density of electrons. As a function of r s 
the ground state of the 2D electron gas can be, for in- 
stance, a Wigner crystal (in dilute limit, at r s > 35 |(J), 
a ferromagnetic Fermi liquid with spontaneous magne- 
tization or a paramagnetic Fermi liquid. An improved 
hypernetted-chain approximation was employed to study 
the 2D electron liquid at full Fermi degeneracy 0,0- The 
results indicate the paramagnetic phase as the ground 
state of the system until the Wigner crystallization den- 
sity, even the separation with polarized states become 
minute (the energy difference with the ferromagnetic 
phase diminishes to milliRydbergs). In contrast, the 
transition from paramagnetic fluid phase to Wigner crys- 
tallization was recently studied for a 2D electron gas by 



diffusion Monte Carlo simulations including the backflow 
corrections [jj showing the stability range of a polarized 
phase in-between. The earlier work including the back- 
flow correlations did not find the polarization transition 
[To). More recently a weakly first-order transition to a 
spin-polarized state was found to occur shortly before 
the Wigner crystallization [TT) . 

Hubbard model is the simplest tool that can capture 
some of the salient properties of correlated systems. A 
delocalized Coulomb phase in 2D was observed in the cal- 
culations of by Waintal et al. 0] which showed the tran- 
sition from the Anderson insulator towards an extended 
phase. Mean-field approximation considered in an unre- 
stricted Hartree-Fock ansatz (UHF) was compared with 
constrained path quantum Monte Carlo method for re- 
pulsive Hubbard model at half- filling 0] and found both 
Anderson and Mott insulator phases. In a 3D system 
numerical mean-field studies by Tusch and Logan 14] at 
the unrestricted Hartree-Fock level revealed the phase 
diagram at half-filling. The phase diagram in a 2D sys- 
tem is studied by Hirsch|l5| using a square lattice with 
nearest-neighbor hopping. 

Nowadays there is also a growing interest in the spin 
polarization and Stoner instability of finite-size disor- 
dered systems O ) E3 ) E3 j 13 > p0|- In quantum dots, 
the singlet-triplet transitions even for weak interactions 
cause switching between states and a kink in the con- 
ductance 0]. For interaction strength below the Stoner 
criterion a nonzero spin ground state was found for an 
ensemble of small metallic grains in the Hartree-Fock ap- 
proximation 0] . Close to the Stoner instability in a 2D 
metallic system there is an exponentially sm all p roba- 
bility for the appearance of local spin droplets 21j. The 
tendency of disorder to exhibit the magnetism of ground 
state was already pointed out in some recent works 18], 
22], 23]. Meanwhile the fluctuations of the interaction 
matrix for different disorder ranges seem to have an im- 
portant role in the spin polarization transition for finite 
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systems. At crossover between random matrix ensembles 
the correlations between different eigenvectors give rise to 
enhanced fluctuations of the interaction matrix elements 
(IME) in small metallic grains and semiconductor quan- 
tum dots 24]. When their fluctuations dominate a 2D 
random interaction model predicts an increased proba- 
bility of a minimal spin for the ground state lowering in 
the tails of the energy band [lTj and for a two-orbital 
model they also help in stabilizing a ground state with 
minimal spin |25j . 

Our goal in this work is to study the influence of dis- 
order and interactions on the polarization of the ground 
state for a 2D finite lattice. We treat the interaction 
effects at the unrestricted Hartree-Fock level (UHF), ne- 
glecting the off-diagonal contribution of the Coulomb ma- 
trix. This approximation permits us to write an effective 
one-body Hamiltonian, suitable for calculating the con- 
ductances in the Landauer-Biittiker formalism [2(|. We 
solve numerically the self-consistent set of equations for 
2D systems with Hubbard interaction grater than the 
Coulomb interaction strength (Ujj S> U). For different 
disorder ranges we study the fluctuations of the inter- 
action matrix elements and point out their role in the 
disorder induced spin polarization. 

The rest of this paper is organized as follows. In the 
next section we outline the theoretical framework of un- 
restricted Hartree-Fock approximation and develop an ef- 
fective Hamiltonian for the problem at hand. In Sec. Ill 
we present our results for the fluctuations in the interac- 
tion matrix elements in various disorder regimes. Section 
IV discusses the possibility of disorder enhanced spin po- 
larization within our model. We close by listing our main 
conclusions in Sec. V. 



II. THE UNRESTRICTED HARTREE-FOCK 
FORMALISM 



We study a disordered rectangular lattice with n 
fermions on L — L x ■ L y sites with vanishing bound- 
ary conditions imposed. If c ia and Ci a are the creation 
and annihilation operators of an electron at site i with 
spin a, the Hamiltonian will be defined as follows: 



H = 1 E c *W + E 



<i,j>,cr 

+u E 



(i) 



where rij )CT = c| Cj ff is the occupation number operator, 
u>i is the random Anderson disorder (ajj S [— W/2, IF/2]), 
UHub and U are the strength of the Hubbard and long- 
range Coulomb interaction, and t is the energy unit of the 
hopping integral (it will be considered equal to unity). 
For uji set to zero, the Hamiltonian of Eq. (1) is the ex- 
tended Hubbard model, used for studying the correlated 
systems |27j. The Hubbard term was first introduced in 
[28j and the argument in favor of it is that only electrons 



with opposite spins can occupy the same state (Pauli ex- 
clusion principle). 

We write the Hamiltonian given in Eq. (1) in a ba- 
sis formed by two sets of orthonormalized single-particle 
states 4>a(i), 4>a(i) with a E [1, L 2 ] (unrestricted Hartree- 
Fock orbitals). Creation and annihilation operators in 
the new basis cJ CT and c Q(T are given by the transformation 

4,<r = Ei<(*K> andc aiCT = Y, l ( l ) a*( i )ciM- Employing 
the variational principle, we look for the self-consistent 
set of eigenvectors <^(i) that minimize the ground state 
energy Eq for a given number of electrons n and a def- 
inite component of spin along an arbitrary z axis, S z . 
The ground state will be a Slater determinant \^q) that 
correspond to the n a electrons in the states </>^(«) for 
a = 1, n a and a =f , \ (rt|+n|=?i): 



(2) 



where |^o) is the vacuum [23]. The eigenvectors 
and the corresponding Lagrange multipliers t G a that give 
the lowest energy solution for the ground state energy 
Eq = ^Iq\H\^g) will be calculated by the variational 
principle. After a straightforward calculation one obtains 
the following spin-separable Hamiltonian H e ^ — H a 
for cr =t,|: 



H ° = E C *W + E N + Vd (0 + V Hub (OK 

<i,j> i 



(3) 



with <j>a(i) and eigenvectors and eigenvalues of Eq. (3): 
Hal^a) — e a\^a)- The time reversal symmetry of H en- 
sures the possibility of choosing a real system of eigen- 
vectors. 

The matrix elements Vo(i), V§ x (i,j) and Vfi ub (i) (di- 
rect, exchange and Hubbard interactions, respectively) 
are related to the single particle eigenvectors by the fol- 
lowing formulae (we stress that in the following a means 
the opposite spin of a): 



Vb(») 



2^E 



V% x (i,j) = 2U 



(4) 



(5) 



Vg ub (i) = U Hub J2\c/>°(i)\ 2 (6) 

a=l 

The ground state energy is given by: 

e g = E e « + C/ E v Zff - u E kB*T 



OLff.B'0 



+ U H ub E C aa,% 



(7) 



a<7,(3<7 



where the second and the third terms are the contribu- 
tions of direct and exchange interaction, respectively, to 
the ground state energy, and the last term is the Hub- 
bard interaction. The summations in Eq. (7) are over 
the occupied states \ao) from Eq. (2). The Coulomb and 
Hubbard interaction matrix elements are given as 



v uo,po' — / j 



c 



E 



(8) 



(9) 



Omitting the off-diagonal matrix elements ((a, (3) ^ 
(7,(5)) in Eq. (7) yields the Hartree-Fock approximation. 
The reason for doing this will be discussed subsequently. 

The basic parameters of our theoretical framework are 
the electron number n, the long-range Coulomb inter- 
action U, the Hubbard interaction strength UHub, and 
the disorder amplitude W. Equations (3-7) constitute 
a set of self-consistent equations to yield the eigenstates 
(j>a , t a a and the energy Eq for a given set of spin occu- 
pation numbers (n|,n^) with n-\ + rij = n. At every 
step of the numerical calculation, the matrix elements 
[described by Eqs. (4-6)] of H e ^ will be calculated for 
the spin occupation numbers that give the minimum en- 
ergy in the previous step. The first Slater determinant of 
Eq. (3) will be formed by noninteracting eigenfunctions 
of Hamiltonian [Eq. (1)] with U = UHub = 0. Varying the 
spin occupation numbers so that M z £ [0, 1] we calculate 
Eq{M z ) whose minimum value gives the spin magnetiza- 
tion M z of the system. The stability of the self-consistent 
solutions is very sensitive at larger Coulomb interaction 
strengths, forcing us to keep a relatively small values of 
the interaction parameter U. Because of the spin rota- 
tional symmetry of the Hamiltonian, one limits the nu- 
merical calculations only to positive spin magnetization 
(n T > nj. 

Increasing the amplitude of disorder the spectrum of 
noninteracting Hamiltonian H = tX}<i j> a c lcr c jc + 
CT Winder evolves from weakly to strongly localized 
regime while the eigenvalue statistics exhibits a crossover 
between Wigner and Poisson distributions j^D] ■ As a con- 
sequence the distribution function of neighboring level 
spacing undergoes a continuous crossover from Wigner 
surmise Pw(s) — tts/2 ■ exp(— irs 2 /4) with the vari- 
ance 5(s) = \/4/ir — 1 ~ 0.52 toward Poisson Pp(s) — 
exp(— s) with 6(s) = 1. (s is the dimensionless level 
spacing measured in the unit of the mean level spacing 
A = (S)). This scenario is also available when the Hub- 
bard interaction is taken into account. For nonzero spin 
magnetization, the whole spectrum of H e ** consists of 
two distinct sequences of the two spectra of the H-\ and 
Hi- The spectrum of H^ depends on the site occupa- 
tion number nj,(i) [c.f. Eq. (6)] that will become zero at 
large Uh (when M z = 1 and nj = n), conserving the 
eigenstates degree of localization as in the noninteract- 
ing case. Meanwhile, the spectrum of H± depends on 
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FIG. 1: The eigenspectrum of H-\ (a) and Hi (b) versus the 
disorder amplitude W for a system L 2 =8-9, Uh = 4, U — 
and the electron number n = 8. The mark symbols between 
the eigenvalue lines show different values of the variance of 
the neighboring level spacing 5(s). From left to right they 
indicate: S(s) < 0.5 (crosses), S(s) £ [0.5,0.54] (times), <5(s) £ 
[0.6,0.7] (stars) and 5(s) > 0.7 (boxes). The averages are 
made over an ensemble of 500 configurations. 



the site occupation number (whose average over 

disorder configurations will become equal to n/L 2 ) and 
moves up at a rate of (V Hub ) disorder = U Hu b ■ n/L 2 . The 
spectrum of H± becomes delocalized as a counterpart of 
V Hub (i) against the disorder potentials w(i). In Fig. 1 we 
show the differences in the two sequences of spectrum, 
Hi states being more delocalized compared to iff states 
at the same W. The variance 8(s) ~ 0.7 in Fig. 1(a) for 
W = 10 as for the noninteracting spectrum in the inset 
of Fig. 2(c). In Fig. 2(c) the variance S(s) > 0.8 for a 
system size L 2 — 8 • 9 with disorder amplitude W = 20. 
One notes that for a finite system the discreteness of the 
energy spectrum does not allow for a real Poisson distri- 
bution and the level repulsion characteristic to Wigner 
surmise is pres ent even at the critical point of Anderson 
transition 36]. The effect of the Coulomb interaction 
(U 0, Uuub = 0) over the spectrum properties con- 
cerns the localization effect of long range interaction for 
£ > L 23] , the edge localization of the occupied Hartree 
levels j3J suggested also in and increasing gap near 
e_F j34j. In the large lattice size limit the Coulomb inter- 
action induces a delocalization process by the meaning of 
Poisson- Wigner transition of level spacing distributions 

m. 



III. INTERACTION MATRIX FLUCTUATIONS 



The Coulomb interaction matrix elements (CME) V2'p 
from Eq. (8) and the Hubbard matrix elements (HME) 
C^'p from Eq. (9) are calculated for disorder functions 
in different localization regimes. To this purpose, we 
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FIG. 2: The mean values of the interaction matrix elements 
versus amplitude of disorder W. L 2 = 8 ■ 8. (a) Direct matrix 
elements V"" with a = 1,2,3 and with a,f3 = 1,2,3; 
(in the inset are represented the same matrix elements for 
L 2 — 14 • 14). (b) Exchange matrix elements with 

a — 1 ... 10 and the off diagonal matrix elements with 
a/3, 7,<5 = 1,2,3. In the inset V£a+2 for the same system, 
(c) Hubbard matrix elements C aa and level spacing S aa +i = 
e a +\ — e Q . In the bottom inset the off diagonal Hubbard 
elements C a /3 with a /3. a,/3 — 1,2,3. In the top inset 
the variance of neighboring level spacing shows the transition 
from extended regime (8(s) = 0.52) to localized one (S(s) — * 
1). The averages are made over 1500 configurations. 
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FIG. 3: The probability density distributions of the inter- 
action matrix elements. The system dimension L 2 =8-8. 
The elements involved and the disorder amplitude are marked 
in every figure. In (e . . . h) we represented the distributions 
of direct and Hubbard interaction matrix elements with any 
a, f3 £ [1, 11]. For W = 10 all distributions (except for the off- 
diagonal element) are asymmetric with the long tails toward 
the values larger then the disorder averages. The histograms 
are represented so that the calculated integral be equal to the 
mean value of the corresponding element. 



use the eigenvectors <p a (a — 1 . . . L 2 ) generated by the 
noninteracting Hamiltonian Hq for different degrees of 
disorder characterized by W. For a 2D lattice with 
L 2 = L x ■ L y = 8 • 9 we have identified three disor- 
dered regimes: (i) extended, when 6(s) ~ 0.52 ± 0.01 for 
W G [1,4] (Fig. 2(c)); (ii) strongly localized regime with 
£ < L and Poisson like repulsion eigenvalues (S(s) > 0.8 
and W > 15); (Hi) weak localization as intermediate 
for W € [4, 15] when the localization length £ crosses 
the system length L. The identification of these three 
regimes of localization is made for the purpose of analyz- 
ing numerical results presented in this work and to help 
in discussing the features of the interaction matrix ele- 
ments values. When transport properties are involved, 
extended and weakly localized regimes have a common 
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meaning, usually addressed as a metallic regime with 
£ > L and g > 0.5 (where g is dimensionless conduc- 
tance). In the following we discuss the main properties 
of the interaction matrix elements in various disordered 
regimes. 

A. The disorder average of the CME with equal in- 
dexes (V^'a) depends on the eigenvectors correlations 

(I0o(*)l * 10a 0)1) w i tn i j- ^ nas large values in 
the extended and weak localization regimes where the 
correlations of the same eigenvectors are proportional to 
1 j g for large distances }33j and decreasing to zero when 
the states became strongly localized on singular lattice 
points (f> a (i) = <5i,i Q (Fig. 2(a)). We note that these ma- 
trix elements together with the Hubbard term C£ give 
the whole interaction energy of the two electrons (with 
opposite spin) in the state a (c.f. Eq. (7)). In contrast, 
the direct matrix elements (V^'p) with a ^ [3 are de- 
creasing when the disorder is increased, keeping constant 
values for strong disorder (Fig. 2(a)). In this case, as- 
suming that the two states </f> a ,0/3 are strongly localized 
in the two random points i a ,i/3 we obtain (V"'^) con fi g — 
(I2i a7 tj0 | j -j | )r.p. i where the last average is made over 
all possible random pairs of points (i a , ip). For a 2D sys- 
tem with I? — 8 • 8 and vanishing boundary conditions 
(Ei^ U^r>r.p. = 0.319 and for L 2 = 14 • 14 it is 
equal to 0.194 and these values are close to the averages 
reached by (V^f) con f ig in Fig. 2(b). 

B. The exchange matrix elements (V^'p) and the off- 
diagonal matrix elements (V^'f) (with (a, /3) ^ 

in Fig. 2(b) have maximum values in the weak localiza- 
tion regime for W ~ 5.5 ± 0.5 when £ > L. This is 
the region where the eigenvector correlations are large. 
For uncorrelated states in the strongly localized regime, 
the off-diagonal matrix elements are decreasing to zero 
faster than the exchange matrix elements. This is also 
possible for weakly localized regime when the exchange 
elements are calculated between the uncorrelated states 
whose eigenvalues are well separated. In the inset of 
Fig. 2(b) the exchange matrix elements are lower 

when eigenvalue differences S a p are increased. The neg- 
ative values in Fig. 2(b) for V£'s when 1 < £ -C L are a 
property of the finiteness of the system and the orthogo- 
nality of the eigenvectors 4> a ,4>p. 

C. The Hubbard matrix elements are given by the in- 
verse participation ratio of a state <fi a , C a>a = Ei^aM 
that evolves from 1/L 2 for extended states to 1 for 
strongly localized ones (Fig. 2(c)). The off-diagonal 
terms C a .p = Ei^aW^lM decay to zero for strongly 
localized states and have maximum values in the weakly 
localized regime where the correlations of two different 
eigenf unctions are proportional to 1/g [33J. This sug- 
gests that the spin polarization transition when Uh i= 
is larger in the strongly disordered regime. As a charac- 
teristic of a finite system, the saturation of the inverse 
participation ratio (C QiQ , = 1) corroborated with the lin- 
ear increase of the mean level spacing makes impossible 




Amplitude of disorder W 

FIG. 4: The spin magnetization (M z ) vs. disorder amplitude 
W for the filling factor n/L 2 = 1/9 and for different inter- 
action strengths (Uh = 4,(7 = 0), (Uh = 2, U = 0) and 
(Uh = 0, U = 0.4). The electron number and the system size 
are n — 4 and L 2 =6-6 (crosses), n = 6 and L 2 =6-9 (times), 
n — 8 and L 2 = 8-9 (stars) and n = 10 and L 2 — 9T0 (boxes). 
In the inset for Uh = 4 and L 2 =8-9 the effect of disorder 
induced spin magnetization disappears when W > 50. The 
averages are made over 200 ~ 1000 configurations. The large 
values of (M z ) for L 2 =6-6 when W < 4 are due to the 
degeneracies of the spectrum and have no significance for the 
disorder induced transition. 

a spin polarization transition for very large disorder. It 
does not take place in the thermodynamic limit (L — > oo 
and lattice parameter a = const.) when the neighboring 
level separation S — > 0. 

The fluctuations of the interaction matrix elements for 
different disorder ranges are presented in Fig. 3. For small 
disorder, the probability density distributions of the IME 
are large and they typically have a Gaussian shape. The 
exchange (and off-diagonal) matrix element fluctuations 
become sharper when their mean values become zero 
and the inverse participation ratio has larger fluctuations 
when disorder is increased, with a shape evolution simi- 
lar to increasing the Coulomb interaction strength (com- 
pare with Fig. 5 in Levit and Orgadj^). We notice the 
asymmetric long tail fluctuations of direct, exchange and 
Hubbard matrix elements (Fig. 3 for W = 10), the last 
two permitting an increased probability of the spin po- 
larization transition even when the Stoner criterion does 
not allow it. 



IV. DISORDER ENHANCED SPIN 
POLARIZATION 

The formula of Eq given in Eq. (7) for the ground 
state energy written as an unrestricted Slater determi- 
nant [Eq. (2)] means that we omit the off-diagonal matrix 
elements from expansion of Eq- In Fig. 2(b) the largest 
values of the off-diagonal elements in the weak localiza- 
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FIG. 5: The spin magnetization (M z ) vs. the filling factor 
n/L 2 . (a) The interaction is Uh =2,U = and disorder W = 
10. In upper inset the transmittance of the last occupied level. 
The three set of lines are calculated for L 2 — 6 • 6 (crosses), 
L 2 = 9- 10 (times) and L 2 = 9- 10 (stars). In down inset (M z ) 
is decreasing to zero at n/L 2 cx 1. (b) For W = 10 and L 2 = 
9 • 10 the interaction strengths are Uh =2,U = 0.1 (crosses), 
U H = 2, U = 0.2 (times) and U H = 3,U = 0.3 (stars). The 
averages are made over 200 ~ 1000 configurations. 



tion regime (V off ~ 0.008 ± 0.01) are much lower than 
the mean level spacing (S(S) ~ 0.2±0.1 in Fig. 2(c)) pro- 
viding us a good approximation in neglecting them when 
the Coulomb interaction strength U<5. The Stoner cri- 
terion for itinerant magnetism is related to the exchange 
energy gain by an antisymmetric spin polarized two-body 
wave function that compensates the reduction of the ki- 
netic energy. With our ansatz, in the regime of extended 
metallic waves it is expected that the large exchange ma- 
trix elements (especially the elements V^_f") can favor 
a singlet-triplet transition when U ■ V ex ~ S. For a 2D 
system with L 2 = 8 • 9, (S) ~ 0.2 and {V ex ) ~ 0.05 (in 
the weak localization regime for W = 6) the interaction 
strength U must be at least 4. Besides, the enhanced 
interaction due to the increased return probability ex- 
hibits a finite temperature partial spin polarization |20| . 



By increasing the disorder, the exchange matrix elements 
have positive long tail fluctuations at values larger than 
the disorder averages (Fig. 3(a)) compensating the level 
repulsion of neighboring levels (and favoring also a less 
probable singlet-triplet transition). This will not be pos- 
sible for strongly localized states when any overlap inte- 
gral is equal to 0. In Fig. 4 we can see this scenario for 
a small value of the Coulomb interaction U = 0.4. Diffi- 
culties in solving the self-consistent equations [Eqs. (3-7)] 
for larger values of the interaction strength U prevent us 
from obtaining reasonable polarization transition when 
Uh = 0. The small values of M z in this case represent 
the rare events in the disorder configuration space when 
the fluctuations of V ex are larger than level repulsion. 

The disorder induced spin polarization was previously 
obtained |22j | for the periodic boundary conditions with 
U = Uh and for a truncated many-body Hilbert space. 
This is possible for nonzero Hubbard interaction strength 
Uh when the evolution of S z is increased by disorder 
(different curves in Fig. 4) . The full polarization regime 
(M z = 1) for a large value of W seems to be stable for 
the constant filling factor n/L 2 when the system size is 
increased. As it was stressed before, for linearly increas- 
ing level spacing at larger disorder values, the spin align- 
ment will be turned off when Uh -C a . a = Uh < S a , a +i ~ 
W/L 2 (see inset of Fig. 4). In Fig. 5 we calculate the evo- 
lution of M z varying the 2D electron density (n/L 2 ) for 
W = 10. When the electron number is increased the sys- 
tem shows a smooth decrease in M z . At the same time, 
the transmittances calculated with Landauer formula for 
electrons at Fermi level show the features of 'insulator- 
metal' transition (in the inset of Fig. 5(a) the transmit- 
tances are enhanced by a factor of about 4) . As a function 
of system properties (W and Uh) the spin magnetization 
can decrease from full or partially spin alignments. Our 
data show a possible transition toward a fully polarized 
state when the electron number n is decreased, sugges- 
tive of and similar to experimental results of Vitkalov et 
at. and Shashkin et al.\^. This takes place only for 
the values of interaction/disorder parameters that allow 
for a total spin alignment at low electron density (see 
Fig. 5(b)). 



V. CONCLUSION 

In this work we have studied the interplay between 
the disorder and interaction effects with regard to the 
spin magnetization of correlated electrons on a finite 2D 
lattice. Our numerical calculations suggest that spin po- 
larization occurs for a finite 2D system with a model of 
disorder both for extended and localized regimes. The 
main results we derive are listed below. 
(i). The level repulsion of neighboring states in the ex- 
tended (and weak localization) regime is balanced by the 
large values of exchange energy or by their long tail asym- 
metric fluctuations when the system is gradually local- 
ized and the eigenfunctions become (extremely) sparse. 
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When the exchange elements decrease towards zero (in- 
creasing W) the spin polarization is not allowed unless 
strong Hubbard repulsion is turned on Uh 7^ 0. 
(ii). The disorder induced spin polarization is calculated 
for a finite 2D strongly correlated system and the data 
are self scaled in the limit of n/L 2 = constant. In the 
framework of unrestricted HF orbitals this is due to the 
increased separation of up and down sequences of the 
spectrum. When the full polarization is assigned, the 
short range Hubbard interaction will have no localiza- 
tion effect over the occupied orbitals while the opposite 
spin orbitals become increasingly delocalized. For a finite 
system the disorder induced spin polarization will be cut 
off when the linearly increasing level spacing exceeds the 
saturation value of inverse participation ratio. 



(Hi). When the number of electrons n decreases the 2D 
system can show a transition to a fully spin-polarized 
state. This, however, takes place only when the disorder 
(W) and interaction (Uh) strengths are properly tuned 
to allow for a ferromagnetic spin alignment in the system 
at low densities. 
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